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Chain-like transitions in Wigner crystals: Sequential or non-sequential? 
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The structural transitions of the ground state of a system of repulsively interacting particles con¬ 
fined in a quasi-one-dimensional channel, and the effect of the interparticle interaction as well as 
the functional form of the confinement potential on those transitions are investigated. Although the 
non-sequential ordering of transitions (non-SOT), i.e. 1-2-4-3-4-5-6-... sequence of chain configura¬ 
tions with increasing density, is widely robust as predicted in a number of theoretical studies, the 
sequential ordering of transitions (SOT), i.e. 1-2-3-4-5-6-... chain, is found as the ground state for 
long-ranged interparticle interaction and hard-wall-like confinement potentials. We found an energy 
barrier between every two different phases around its transition point, which plays an important 
role in the preference of the system to follow either a SOT or a non-SOT. However, that preferential 
transition requires also the stability of the phases during the transition. Additionally, we analyze the 
effect of a small structural disorder on the transition between the two phases around its transition 
point. Our results show that a small deformation of the triangular structure, change dramatically 
the picture of the transition between two phases, removing in a considerable region the non-SOT 
in the system. This feature could explain the fact that the non-SOT is, up to now, not observed in 
experimental systems, and suggests a more advanced experimental set-up to detect the non-SOT. 

PACS numbers: 81.30.-t, 37.10.Ty, 82.70.Dd, 52.27.Lw 


I. INTRODUCTION 

At low temperatures, a classical system of charged 
particles arranges itself in a close packet structure, also 
known as Wigner crystalid— . This organization of parti¬ 
cles allows the lowest energy configuration and, due to 
the absence of kinetic energy, the arrangement of parti¬ 
cles results in a stable crystalline structure. In a two- 
dimensional system, Wigner crystals have a hexagonal 
lattice structurei^^— . 

For a quasi-one-dimensional (QID) system of classical 
charged particles confined in a parabolic channel, Pia- 
cente et alX predicted a non-sequential ordering of tran¬ 
sitions (non-SOT) between ground state (GS) configura¬ 
tions: 1-2-4-3-4-5-6-... chain-like structures with increas¬ 
ing particle density. They found that the range of the 
interaction between particles does not affect the ordering 
of the transition, when the particles are confined by a 
parabolic potential^^— . The structural transition from 
two- to four-chain configuration (2 —>■ 4) occurs, in the 
case of a non-SOT, through a zigzag transition of each of 
the two chains and a simultaneous small shift along the 
chain, which makes it a discontinuous transition^. 

Although this non-SOT has been found as the GS in 
a number of theoretical work s it was not, to the 
best of our knowledge, observed in experiments so far. 
Instead, in the experiments they observed a direct tran¬ 
sition from two- to three-chain configuration (2 —>■ 3), 
allowing the system the follow a usual sequential or¬ 
der of transitions (SOT), with increasing linear den¬ 
sity, as reported in the case of electrons on liquid He 
at low temperaturesi^^— and even in dusty plasma clus¬ 
ters!^. The SOT in a system of electrons on liquid He 
has been indirectly measured by interpreting the step¬ 
like increment of the conductance of a channel, while 
an electrical force guides the motion for the particles 
in the structur o^^d^ . Theoretical works have modeled 


this system, evidencing the SOT during that dynamical 
processii^— , although modifying the shape (i.e., increas¬ 
ing the length) of the confining constriction was shown to 
facilitate the observations of the non-SOT—. Similarly, 
a SOT has been predicted theoretically for a binary mix¬ 
ture of repulsive particles trapped in a channel^^, and 
for an Abrikosov-vortex arrangement in a superconduct¬ 
ing slab for low temperatures^i^— , as well as for Pearl 

vortices^i. 

Ikegami et ah^ observed a periodic change of melted 
and ordered states as a function of linear density, i.e. 
re-entrant melting. However, in that work melting tem¬ 
perature was just a measure of disorder of the system. 
In terms of structural transitions, they always observed 
SOT, in agreement with Ref. Q if temperature is above 
a certain value. However, following the results of Pia- 
cente et al^ one might expect that lowering temperature 
would result in the non-SOT, which was not seen exper¬ 
imentally within the attainable temperature rang o^^d^ . 
Thus the important question rises: how to optimize ex¬ 
perimental set-ups in order to realize the observation of 
the non-SOT? 

Previous theoretical works have shown the robustness 
of the non-SOT for a system of particles trapped in a 
perfect parabolic channel irrespective of the range of the 
interaction between the particles^, and also for a system 
of particles with a fixed interparticle interaction, confined 
in QID channels with potentials of various functional 
formj^i However, it has been recently proven that it is 
possible to hnd a SOT region by tuning the parameters 
of the confinement potential, e.g. by invoking a profile 
similar to the Bean-Livingston barrier for vortices in a 
superconductor—. At first glance, we note that two ex¬ 
treme cases of this problem have been analyzed: a fixed 
confinement varying the range of the interaction, and a 
fixed interaction range varying the profile of the confine¬ 
ment channel. From these results a trend of the ordering 
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of the transitions was revealed, but it still did not resolve 
two important issues: why the system prefers following a 
non-SOT rather than SOT? And, most importantly, why 
experiments did not find a non-SOT if it is a preferable 
scenario? In this work, we present an in-depth analysis 
which answers these questions. 

The present paper is organized as follows. We first 
give, in Sec. m an overview of the model and the per¬ 
spective of the current work done in this area. Sec. uni 
is devoted to the analytical study of the effect of the 
interaction range on the ordering of the structural tran¬ 
sitions, for different confinement potentials. In Sec. IIVI 
we induce transitions from two- to three-chain (2 —>■ 3) 
and from two- to four-chain (2 —> 4), and calculate the 
energy barrier for each transition, analyzing the stability 
of the configurations formed in the vicinity of the tran¬ 
sition point. The effect of small imperfections on the 
energy barrier and on the ordering of the transitions are 
analyzed in Sec. [V] Finally, our conclusions are given in 
Sec. EH 


II. MODEL SYSTEM 


First, we consider a two-dimensional infinite system 
of identical interacting particles with electric charge q 
and mass m, which move in the xy-plane. The parti¬ 
cles are confined by a one-dimensional potential limit¬ 
ing their motion in the y-direction, forming a quasi-one- 
dimensional channel along the a;-axis [Vconfiu)]- The to¬ 
tal energy of the system is given by: 

OO OO OO 

H — Vint{\^i — Tj I) -f VconfiUi), (1) 
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where is the relative position of the i-th particle in 
the system, while Vintir) represents the pairwise inter¬ 
particle interaction, which is taken as a screened power- 
law potential, which will allow for the simulations of both 
short- and long-range interactions, as follows: 
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where the parameters A and n allows us to tune the range 
of the interaction between particles in the system, e is 
the dielectric constant of the medium the particles are 
moving in, and R is an arbitrary length parameter which 
we introduced to guarantee the right units. 

Due to the importance of the profile of the confinement 
potential on the ordering of the transitions, as shown in 
Ref. [13, we use the following two different functional 
form of the confinement potential, which allows us to vary 
the profile of the channel continuously from a parabolic¬ 
like to a hard-wall potential: 
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where /3 and 7 control the sharpness, v the strength, and 
2/0 the effective width of the confining channel. In the 
following, we refer to Va(/ 3,2/) [44(7,2/)] as exponential 
[Gaussian] confinement. 

In dimensionless form, the interaction and the confine¬ 
ment potentials in our model become: 


5 


(5) 

cosh(/3?/) — 1 
cosh(/3yo) - 1 ’ 


( 6 ) 

g-7^(y-yo)^ _l_ g-7^(y+yof 

5 

(7) 

is expressed in units 

of 


Vint (r) — " 

VA{l3,y) = v^yl 
VBh,y) = (r'^yl 


where the energy 

Eq = (mwg/2)"/4+2)(g2/g)2/(n-H2)_^2(n-l)/(n+2) 

and all distances are expressed in units of 
ro = (2gVma;ge)i/4+2)^(n-i)/(n-H2)^ dimen¬ 

sionless frequencies are given by t> = va/ujq and 
O’ = vbI^Qi while ojQ measures the strength of the 
confinement potential, and the screening of the pairwise 
interaction is k = ro/A. The dimensionless linear density 
77 is defined as the number of particles per unit of length 
along the unconfined direction. 

Concerning the sequence of the GS configuration for 
increasing system density, previous works have shown 
that: 1) In case of parabolic confinement, the GS follows 
the non-SOT irrespective of the range of the interparti¬ 
cle inter action!. 2) In case of the interaction potential 
given by n = 1 and k = 1 in Eq. ®, and varying the 
profile of the channel, the GS follows the non-SOT in the 
limiting cases of hard-wall and parabolic-like profile^!. 
The latter study also showed that non-SOT is present in 
all systems for intermediate values of the shape param¬ 
eter, evidencing that the non-SOT is extremely robust 
for a broad range of possible profiles and shape parame¬ 
ters. However, in case of Gaussian confinement, SOT for 
the GS was found within a small window of the shape 
parameters!!, indicating that the shape of the channel 
profile is of crucial importance when one is looking for a 
SOT in the system. 

These results bring in evidence the strong relation be¬ 
tween the range of the interparticle interaction and the 
confinement profile. Due to the fact that in experiments 
it is very common to observe the SOT but not the non- 
SOT, a complete and detailed study about the effects of 
the range of the interparticle interaction and the profile 
of the confinement channel on the ordering of the tran¬ 
sitions, are needed, as well as a trustable recipe for the 
experimentalist about the regions where the non-SOT 
could be observed. The present paper will address those 
problems. 


III. INTERACTION RANGE VS. 
CONFINEMENT PROFILE 


Analytical calculation of the energy of a system of dif¬ 
ferent chain-like structures, is performed by following the 
model proposed in Ref. [ 23 . As a result, we found ’’shape 


2 











3 


c 

o 

o 

*■ 4 —» 

c 

o 

c. 

o 

Q. 

X 

lU 


c 

o 

o 

c 

CC 

'(/) 

V) 

CC 

CD 



0 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 

ri ri ri 


FIG. 1. (Color online) Phase diagrams of the ground state for exponential (upper pannel) and Gaussian confinement (lower 
panel), for different interaction potentials as indicated in the figure. Each phase diagram shows the GS as function of shape 
parameter (/3 or 7 ) and linear density of the system {rj). Solid and dashed lines represent first and second order transitions, 
respectively. Regions highlighted by a line pattern indicate the existence of SOT in the GS of the system. 


parameter vs. density” phase diagrams for different con¬ 
finement potentials and different ranges of the interac¬ 
tion between particles, as shown in Fig. [T] The phase 
diagrams are the zero temperature GS configuration as 
function of the shape parameters and the linear particle 
density rj. All the transitions are of first order, except 
the zigzag transition between one- and two-chain which 
is of second order. Interestingly, these results show the 
appearance of a region where the GS is guided by a SOT, 
for three of the model systems, as highlighted by the pat¬ 
terned region in each phase diagram. These results evi¬ 
dence the complementing behavior between the confine¬ 
ment and the interaction potential. Previouslj^, SOT 
had been reported theoretically only in case of Gaussian 
confinement for the interaction parameters n = 1 and 
K = 1. In that case, the SOT was found within a small 
window of 7 (including 7 = 1 which models the Bean- 
Livingston barrier for superconducting vortices). 

An extended phase diagram for the ordering of the 
transitions is plotted in Fig. [21 as a function of the pa¬ 
rameters of the interparticle interaction k and n, which 
control the range of the interaction (see Eq. ([5])). The 
solid symbols (triangles for exponential and circles for 
Gaussian confinement) indicate the case when the GS 
configuration always follows a non-SOT irrespective of 
the shape parameters, while the open symbols indicate 
the case when the SOT has been found. From Fig.|2]one 
can see that, the non-SOT is robust for a vast range of 
shape parameters of both confinement potentials stud¬ 
ied. This behavior allows one to interpret the non-SOT 
as the ’’natural” mechanism, which determines the be¬ 


havior of a system of classical particles with increasing 
particle density, irrespective of the interaction between 
them or the shape of the confinement potential of the 
channel. This result is a generalization of previous theo¬ 
retical work a^d°i^'^ to a large set of different interparticle 
interactions and confinement potentials. 



FIG. 2. (Golor online) Extended phase diagram of the ground 
state transition for exponential (blue triangles) and Gaussian 
confinement (red circles), as a function of the interaction pa¬ 
rameters n and k. The solid symbols indicate that GS always 
follows a non-SOT by increasing density, while open symbols 
indicate the cases where a SOT region has been found. 
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(a) Induced transition from two- to three-chain 
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(b) Induced transition from two- to four-chain 
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FIG. 3. (Color online) Paths of the motion for the dragged 
particles, during the induced transition: a) 2 —>■ 3, and b) 
2 ^ 4 . 


Note that the SOT has been found in some regions of 
the shape parameters, when the interaction potential is 
long-range. However, the range of the interaction to ob¬ 
serve a SOT must be even longer in case of exponential 
confinement as compared to Gaussian confinement. This 
finding is very interesting because the experimental ob¬ 
servation of SOT in systems of electrons floating on the 
surface of liquid He, assume a Coulomb-like interaction 
between particles^^, which is clearly a long-range inter¬ 
action. 


IV. ENERGY BARRIER AT THE TRANSITION 
POINT 

In order to understand why the system prefers to follow 
one ordering (SOT versus non-SOT) instead of the other 
one, we compare the energy barrier the system has to 
overcome during the transitions 2 —3 and 2 —4. 

For this purpose, we will focus on the case of the in¬ 
teraction parameters n = 2 and /c = 1, for the Gaussian 
confinement potential given by Eq. because in this 
case the SOT is present in a small window of values of 
the parameter 7 . In order to calculate the energy barrier 
between the configurations, we set the two-chain struc¬ 
ture, which was found as metastable or GS configuration 
close to the transition point (? 7 t), as initial configuration. 
We induce the transitions 2 —3 and 2 —4, at a given 
density {ijo < rjt), by forcing the displacement of some 
selected particles, allowing the rest of the particles to 
adjust themselves in order to reach the most energeti¬ 
cally favorable configuration. The selected particles are 
dragged along straight trajectories as shown in Fig. [SKa) 
and Fig. |3l(b) for the transitions 2 —>■ 3 and 2 —>■ 4, corre¬ 
spondingly. The straight paths are expected to be close 
to the real trajectories of the selected particles taking 
into account the symmetry of the initial and the final 
configurations during the transitions 2 —3 and 2 —>■ 4. 


Following the displacements shown in Fig. |3l we found 
that an energy barrier is formed very close to the transi¬ 
tion point (? 7 o .$7*), and at rjt, the barrier for transition 
2 —^ 3 is always higher than the barrier for transition 
2 —4, as shown in Fig. |4l This behavior shows that, 
in terms of energy cost, the system always prefers the 
transition 2 —>■ 4 rather than the transition 2 —)• 3. This 
result provides an understanding, why the system prefers 
to follow the non-SOT which appears to be a straightfor¬ 
ward mechanism for the transition from a two-chain to 
a four-chain configuration, due to the minimum energy 
cost for the system. 

Now, we extend the calculation of the induced transi¬ 
tion not only around rjt, but to the whole region of pa¬ 
rameters. For every displacement during these induced 
transitions, we calculate the dynamical matrix of the sys¬ 
tem, and by using the Newton optimization^, we calcu¬ 
late the most energetically favorable configurations, to¬ 
gether with their vibrational eigenfrequencies. From the 
symmetry of our system, we know that it has one trans¬ 
lational symmetry, which is numerically evidenced by the 
existence of only one vanishing eigenfrequency. This con¬ 
dition is important in our case, due to the fact that it 
allows us to define the stability of the different configu¬ 
rations found, after the induced transition is performed. 
In Fig. [ 5 ] we show, with symbols, the regions where the 
induced transitions 2 —>■ 3 and 2 —4 are stable; at the 
same time we plot the analytical results for the phase 
diagram of the GS (see Fig. [T]). 

From Fig. [S]one can see that, while the induced transi¬ 
tion 2 —>■ 3 is stable irrespective of the value of 7 , the sta¬ 
bility of the transition 2 —>• 4 is restricted to two different 
regions of 7 . We can observe that our numerical results 
show a perfect match with the transition point analyti¬ 
cally calculated in Section Hill It is worth noting that the 
SOT formed in the window (0.7 ^ 7 1), as found an- 
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FIG. 4. (Golor online) The energy barrier for the induced 
transitions as a function of the displacement of the dragged 
particles, for different values of the confinement parameter 
and electron density. Red dashed line represents the tran¬ 
sition 2 —>■ 3, while blue solid line represents the transition 
2-^4. 
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FIG. 5. (Color online) Stability region of the transitions 2 —>■ 4 
(solid circles) and 2 —>■ 3 (open squares). The regions defined 
by the different colors, represents the analytically calculated 
phase diagram for the GS, as a function of the shape param¬ 
eter ( 7 ) and the linear density ( 77 ), as shown in Fig. [T] We 
took the parameters n = 2 and ac = 1 . 


alytically, arises from the fact that the transition 2 —>■ 4 
is not stable at the intermediate points (7 = 0.8 and 
7 = 0.9). Summarizing, the system ’’naturally” evolves 
following the non-SOT with increasing density, since, in 
all the cases the transition 2 4 is more energetically fa¬ 

vorable than the transition 2—^3. However, the window 
of non-stability of the configurations during the transi¬ 
tion 2 —4, opens the door for the transition 2 —3 as 
the GS, and therefore the SOT takes place. 


V. THE ROLE OF FLUCTUATIONS 

In the above analysis, we assumed that there is no 
disorder. As long as a small amount of disorder does 
not affect the above findings, these can be considered as 
robust and reliable. The goal of this section is to analyze 
the effect of disorder. 

Typically, disorder is produced by thermal fluctua¬ 
tions or by imperfections of the geometry of the channels. 
While the geometry can be made nearly perfect, thermal 
fluctuations are inevitable. Also, when thinking of exper¬ 
imental measurements of electrons moving on the surface 
of He in micro-channels, we should keep in mind that the 
experimental methodsi^ do not allow to control the elec¬ 
tron structure itself. The structure is detected indirectly, 
via measuring the electron current through a narrow con¬ 
striction. Therefore, one cannot judge whether the elec¬ 
trons were perfectly ordered or not before they entered 
the constriction. Of course, the limiting cases of a crys¬ 
tal and liquid can be distinguished by the appearance of 
typical oscillations in the /U-curves in case of a Wigner 
crystali^^— On the other hand, one can expect that 
weak disorder cannot be detected in the electron current 
measurements, i.e., these measurements cannot distin- 


(a) 




Otdy 


9 -•- 9 9 

9 9 9 


•tb 

otdy 


“i; 


• • 


FIG. 6 . (Color online) Schematic representation of the struc¬ 
tural imperfections induced in the structures of (a) three- and 
(b) four-chain configuration. 


guish a perfect Wigner crystal from a slightly disordered 
one. However, it is not known whether the order of tran¬ 
sitions studied above is sensitive to weak disorder. 

In order to analyze the effect of disorder, we intro¬ 
duce a small displacement of one particle (per simula¬ 
tion cell) from its symmetric GS configuration, as shown 
in Fig. |6l Note that in the general case this procedure 
is equivalent to the effect of non-zero temperature (see, 
e.g., Ref. 0 )- In particular, a displacement of 10% of the 
lattice constant from its equilibrium position, according 
to the Lindemann criterion, is treated as (local) melting 
of the crystal. 

Following the procedure described above, we increase 
the particle’s density and analyze the structural transi¬ 
tions. We calculate the relative energy barrier of the tran¬ 
sition from two- to the imperfect three-chain (2 —3*) 
and to the imperfect four-chain (2 ^4*) configurations. 
Fig. [7Ka) shows this barrier plotted as a function of 7 at 
the transition point, circles and squares represent the en¬ 
ergy barrier for 2 —>■ 3* and 2^-4*, respectively. One can 
see that the energy barrier of 2 3* decreases as com¬ 

pared to the one of 2 —>■ 4*, resulting in the shrinking of 
the region where the transition to three-chain is more fa¬ 
vorable than to the four-chain configuration. Therefore, 
the GS of the system follows the SOT in a wide range 
of parameters, as clearly shown in Figs. [7Ka,b). In the 
highlighted gray region in Fig. [TKa), the SOT guides the 
GS due to the instability of the transition 2—^4 config¬ 
uration, as discussed in previous section. However, just 
a small imperfection, as the one modeled in the present 
analysis, breaks the robustness of the transition 2 —s- 4, 
allowing the system to evolve following the SOT for a 
region of 7 much wider than in case of perfect chain-like 
structures (cf. Fig.^]). 

Using these results we re-plot the phase diagram for 
the above studied system (n = 2 and k = 1 ) i.e., shown 
in Fig.[5l In the new phase diagram (Fig. [7Kb)) the GS of 
the system follows the SOT for a broad range of densities, 
starting from 7 Ri 0.7 and higher. Therefore, even a small 
amount of disorder changes dramatically the picture of 
the ordering of the transition in the GS configuration 
of the system. The non-SOT earlier considered as a “ro¬ 
bust” feature (since it survived under various transforma¬ 
tions of the confinement potential and the inter-particle 
interaction) turns out to be fragile against a small de¬ 
formation in the lattice, for low particle density. This 
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FIG. 7. (Color online) (a) Relative energy barrier as a func¬ 
tion of the shape parameter ( 7 ) at the transition point. Open 
circles (squares) represent the energy barrier of transition 
2 —>■ 3* (2 —>■ 4*). The gray highlighted region indicates the 
values of 7 where the transition 2 ^ 4 is not stable, as shown 
in Fig. [5] (b) The same phase diagram as in Fig. [5] but with 
the effect of the induced imperfection in the lattice. 


fragility of the non-SOT explains, why experimentally 
the non-SOT has never been observed under real exper¬ 
imental conditions (e.g., non-zero temperature). 


VI. CONCLUSIONS 

The existence of the so-called non-sequential ordering 
of structural transitions (non-SOT), when increasing the 
particles density, in a system of repulsive particles con¬ 
fined in a quasi-one-dimensional potential was theoret¬ 
ically predicted in a number of works. The non-SOT, 
when the number of chains follows the sequence 1-2-4- 
3-4..., was shown to be the “natural” sequence of tran¬ 
sitions for a broad range of the interaction and confine¬ 
ment parameters. However, in spite of the theoretically 
predicted robustness of the non-SOT, experiments with 
various interacting particles including electrons on the 
surface of liquid He in micro-channels, colloids in narrow 


channels, and superconducting vortices in narrow stripes, 
did not reveal the non-SOT. Instead, the transitions fol¬ 
lowed the SOT, i.e., 1-2-3-4... 

In an attempt to investigate this controversial behav¬ 
ior, we studied in detail the effect of boundaries on the 
sequence of the structural transitions^!. It was shown, in 
particular, that in case of a Bean-Livingston-type barrier, 
which appears in case of vortices in a superconductor, 
one observes SOT rather than non-SOT. In this work, 
we took a deeper insight in this problem, by generalizing 
our study to interparticle interactions of very different 
range. In addition, we analyzed the effect of fluctuations 
which provided us with a deeper understanding why real 
systems display SOT rather than the theoretically pre¬ 
dicted non-SOT. 

In particular, we investigated the structural transitions 
for the ground state (GS) of a classical system of parti¬ 
cles confined in a channel. The inter-particle interaction 
was modeled as a screened power-law potential, and the 
prohle of the channel confining potential was modeled by 
two different functional forms, both being gradually tun¬ 
able from a parabolic- to a hard-wall-like confinement. In 
this context all the GS configurations are presented by 
chain-like structures, which correspond to Wigner crystal 
structures. 

The effect of the interaction between particles on the 
ordering of transitions of the GS configurations, was in¬ 
vestigated for different values of the interaction param¬ 
eters, varying from a short-range to a long-range inter¬ 
action potential. We analytically calculated the energy 
of the chain-like structures at zero temperature, and an¬ 
alyzed the results obtained for different confinement po¬ 
tentials. We found that the non-SOT is present in all 
cases studied and it was found irrespective of the range 
of the interaction and for different energy profiles of the 
confining channel. 

On the other hand, our calculations show that a region 
where the SOT guides the transitions of the GS by in¬ 
creasing particle density, emerges in case of a long-range 
interaction potential (i.e., when approaching the limit 
of unscreened Coulomb interaction), and this region is 
extended for longer-range interactions: The longer the 
range of the interaction, the higher probability to find 
the SOT as a sequence of the GS of the system. 

Two different types of trapping potentials, which con¬ 
trol the profile of the channel where the particles are 
confined, were studied: exponential and Gaussian con¬ 
finement. In both cases the shape parameter allows us 
to control the shape of the channel, varying it from a 
parabolic-like to a hard-wall-like channel by increasing 
the shape parameter. We found that for a hard-wall-like 
confinement, the transitions between different configura¬ 
tions occur at lower particle linear density, and the or¬ 
dering of the transitions is related to the range of the 
interaction potential. Oppositely, in case of a parabolic¬ 
like confinement, we found in all studied cases, that the 
transitions between the different phases are always given 
by a non-SOT. 

The existence of the sequential or non-sequential order- 
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ing of transitions between phases in the GS is determined 
by the existence or absence of the four-chain GS config¬ 
uration “inserted” in the direct transition 2 —> 3. We 
investigated the ability of the system to follow either the 
2 —>■ 3 transition or the 2 —4, by calculating the en¬ 
ergy barrier which the system has to overcome in order 
to reach the Hnal state. We found that, irrespective of 
the confinement profile and the interaction potential, the 
barrier for the 2 —>■ 3 transition is higher than that for the 
2 —>■ 4, thus making the non-SOT as the most “natural” 
ordering of transitions for the GS of the system. 

However, when analyzing the stability of the configu¬ 
ration during this transition (2 —>■ 3), we found that if 
this transition is not allowed as the GS, the reason for 
this behavior is the instability of that state. Our numer¬ 
ical results obtained from the analysis of the stability of 
the states perfectly confirm the behavior of the transi¬ 
tions and the phase diagram calculated analytically for 
this system. The analysis of the stability explains that 
the robustness of the non-SOT is not just due to the ener¬ 
getically favorable arrangement of the particles, but due 
to the instability of the transition 2 —3. 

As further follows from our analysis, the effect of weak 
fluctuations in the chain-like structure can lead to the 
SOT behavior. Thus our simulations showed that even a 
small imperfection could change dramatically the order¬ 
ing of transitions facilitating the appearance of the SOT 
and considerably increasing the window of parameters to 
observe it. This results in a remarkable finding: The ro¬ 


bustness of the non-SOT in the system can be eliminated 
by a weak disorder of the Wigner crystal, e.g., due to a 
non-zero temperature. 

Our results manifested that small fluctuations destroy 
the intermediate four-chain configuration (i.e., the hall¬ 
mark of the non-SOT), facilitating the experimental ob¬ 
servation of the SOT instead of the earlier predicted non- 
SOT as the most “natural” sequence of transitions. At 
the same time, an important result of our analysis is that 
we found the window of parameters where the non-SOT 
still can be potentially found, even in the presence of 
weak fluctuations (e.g., for non-zero but relatively low 
temperature). In particular, our analysis showed that 
the non-SOT is “protected” against small fluctuations in 
the regime of high linear density and smooth confine¬ 
ment, therefore, we expect that the non-SOT would be 
experimentally observed in that region, rather than for 
low density of particles confined in a hard-wall channel. 
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